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^ ABSTRACT 

O 

j—j Context. Rigorous modeling of the Ban /14554 formation is potentially interesting since this strongly po- 
larized line forms in the solar chromosphere where the magnetic field is rather poorly known. 
Aims. To investigate the role of isotropic colUsions with neutral hydrogen in the formation of the polarized 

" ^ Ban /14554 line and, thus, in the determination of the magnetic field. 

00 Methods. Multipole relaxation and transfer rates of the d and p-states of Ban by isotropic collisions with 
neutral hydrogen are calculated. We consider a plane parallel layer of Ba n situated at the low chromosphere 
and anisotropically illuminated from below which produces linear polarization in the /14554 line by scatter- 
ing processes. To compute that polarization, we solve the statistical equilibrium equations for Ban levels 

OO including colUsions, radiation and magnetic field effects. 

o 

• • Results. Variation laws of the relaxation and transfer rates with hydrogen number density nn and temper- 



ature are deduced. The polarization of the /14554 fine is clearly affected due to isotropic colUsions with 
neutral hydrogen although the collisional depolarization of its upper level ^P3/2 is negligible. This is be- 
cause the aUgnment of the metastable levels ^£>3/2 and ^D}/2 of the Ban are vulnerable to colUsions. At 
the height of formation of the /14554 line where mh ~ 2 x 10''* cm"^, we find that the neglecting of the 
coUisions induces inaccuracy of ~ 25% on the calculation of the polarization and ~ 35 % inaccuracy on 
microturbulent magnetic field determination. 

Conclusions. The polarization of the /t4554 line decreases due to collisions with hydrogen atoms. In addi- 
tion, during scattering processes collisions could change the frequency of the Ban photons. To quantitatively 
study this line, one should confront the problem of development of general theory treating partial redistri- 
bution of frequencies and including transfer and relaxation rates by collisions for a multilevel atom with 
hyperfine structure. 
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1. Introduction 

Depolarization due to the Hanle effect can be used to retrieve information on solar magnetic 
fields (e.g. Trujillo Bueno 2001). This magnetic depolarization is usually mixed with a colli- 
sional depolarization. To quantitatively use the Hanle effect as a technique of investigation of 
solar magnetic fields, one needs to know the coUisional rates. In addition, one should include 
properly both the collisional depolarizing and transfer rates in the line formation modeling. For 
instance, Derouich et al. (2007) found that the effect of the collisions, in typical solar conditions 
of temperature and hydrogen density, is particularly important for the very important Ti i /14536 
line- a collisional depolarization of more than 25 % is obtained. In fact, at a first glance, because 
the inverse lifetime of the upper level is larger than the value of the elastic depolarizing rate 
D^, one might think (wrongly) that the effect of the collisions is negligible for the Ti i /14536 
line. Collisional depolarization of this line is mainly due to transfer rates which are generally 
neglected (see Eq. 14 of Derouich et al. 2007). 

Since their ionization potential is rather low, most Barium atoms are ionized throughout the 
low chromosphere (Tandberg-Hanssen & Smythe 1970). The core of the Ba ii /14554 line is 
formed at about 800 km above the photosphere (e.g. Uitenbroek & Bruls 1992). The wings are 
formed in deeper layers of the photosphere. Recently, a renewed interest in the measurement 
and the physical interpretation of the linear polarization of the Ba ii /14554 line have emerged 
(e.g. Malherbe et al., 2007; Belluzzi et al. 2007; Lopez Ariste et al. 2008). It has been shown 
that during the formation of this line the hyperfine structure and the partial redistribution of 
the frequencies effects are very important. In this work, we investigate fully the importance of 
isotropic collisions with neutral hydrogen in its modeling (Sects. 2 and 3). Sect. 4 is dedicated to 
the calculation of the impact of neglecting collisional effects on the magnetic field determination. 
Our concluding remarks are given in Sect. 5. 

2. Modeling of the creation of the polarization in the Ba ii levels 

Let us consider a plane-parallel layer in the low chromosphere (~ 800 km above the photo- 
sphere) formed by Ba n ions and illuminated anisotropically by the photospheric continuum 
radiation field. In a weakly polarizing medium like the solar atmosphere one can safely neglect 
the contribution of polarization to the excitation of Ba n ions. Assuming that the incident radia- 
tion has cylindrical symmetry aroimd the local solar vertical through the scattering center, only 
the multipole orders k = Q (mean intensity J^) and k = 2 (radiation tensor associated to the 
anisotropy) are needed to fully describe the incident radiation. For each wavelength, the value of 
the anisotropy factor (w = "V^/q/j^) and of the number of photons per mode (n = J^{c^ /2hv^)) 
are obtained from Fig. 2 of Manso Sainz & Landi DeglTrmocenti (2002). In particular, for the Ba 
II /14554 line, we obtain w = 0.16 and h = 2.59 x 10"^ which are similar to the values calculated 

* Also associated researclier at CNRS UMR 8 112 - LERMA, Observatoire de Paris, Section de Meudon, 
F-92195 Meudon, France 



M. Derouich: Scattering polarization of the Ban /14554 line 3 

by Belluzzi et al. (2007). We verified that a reasonable modification of w and h (i.e. by around 
10%) do not affect the results of the present work. Obviously, a realistic simulation of the line 
formation conditions requires a careful consideration of the transfer of radiation in a medium that 
is not optically thin. This problem is out of the scope of this paper. The formulae we use to com- 
pute the linear polarization for the case of a tangential observation in a plane-paraUel atmosphere 
is presented, for example, in Sect. 2 of Derouich et al. (2007). 

A simplified atomic model approximation (Fig.[l]) can reduce considerably the numerical and 
theoretical calculations especially if the transfer of radiation effects are considered. However, the 
hypothesis of neglecting metastable t/-sates in the Ba n /14554 line modeling is highly question- 
able. We calculate the polarization p - Qjl using the simplified atomic model (see Fig. [l]) and 
the 5 levels-5 lines model (see Fig.|2| which accounts for the metastable t/-levels. We find that 
the simplified model overestimates the polarization degree of the Ba n /14554: 

^P^5 levels-5 lines model ^ q y (^ y) 

(^'-'simplified model 

Thus, a more realistic diagnostic of the polarization of this line should be performed in the frame- 
work of 5 levels-5 lines model. In the description of the Ba n, we neglect the contribution of the 
hyperfine structure. There are no additional conceptual difficulties for including the contribution 
of the hyperfine relaxation and polarization transfer rates by isotropic collision. Our main con- 
clusions concerning the evidence for collisional depolarization of the Ban /i4554 line remain 
unchanged. 



3. Contribution of collisions in the framework of the 5 levels-5 lines model 

3.1. Definitions of the depolarization, relaxation, and polarization transfer rates and 
statistical equilibrium equations 

The colUsional evolution of the density matrix components Pg(J) is due to the gain-terms denoted 
as polarization transfer rates and to the loss-terms denoted as relaxation rates. The polarization 
transfer rates, denoted by Cj(k J ^ k J', T) if Ej > E'j and Cs{k J <— k J') if Ej < E'j, 
are due to inelastic and super-elastic collisions respectively. E j being the energy of the /-level. 
The left arrow (<— ) is used because these rates represent gain-terms in the statistical equilibrium 
equations (SEE) describing the evolution of the /-level. 

The relaxation rates are the sum of the term Z)*(7, T) exclusively responsible for the de- 
polarization of the level (J) and the terms Cj{k - QJ^k-QJ', T) (Ej < E'j) and 
Cs(k = J ^ k = J', T) {Ej > E'j) corresponding to the population transfer between 
the levels J and J' . We notice that, in principle, the denotations "inelastic" and "superelastic" 
are used for transitions between different electronic sates. Within the same spectral term (inside 
a given electronic state), collisions are called "elastic" even these involving transitions between 
two different /-levels. However, we adopt the denotations Cs and C/ since they are usually used 
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Fig. 1. Partial Grotrian diagram of Ba ii showing the levels and the spectral wavelengths taken 
into account in the case of the simplified model. Note that the level spacings are not to scale. 



in astrophysic^ We denote by T) the relaxation rates of rank k given by: 

R\J, T) = D^{J, T) + R°(J, T) = dHJ, T) + ^ ^2jt \ Cs(k ^ J ^ k ^ Ji, T) (2) 



The overall population of the level J is not affected by purely elastic collisions, r)=0. The 

SEE are: 

[| Pg(J)i^ii = -R^'^J^ T) p^iJ) + ^ C,(k J^kJi, T) p\{Ji) 

h 

+ Y,Cs{kJ^kJ,,T)p\{J,) (3) 

Despite -k < q < k, due to the isotropy of the collisions with neutral hydrogen, the relaxation 
and polarization transfer rates are calculated only for ^ = since they are ^'-independent. The 
effect of the isotropic collisions is the same for density matrix elements with ^ = and for 
off-diagonal ones corresponding to coherences where q + Q. 



' Apart a multiplicity factor of we use the same definitions of the monograph by Landi 

Degl'Innocenti & Landolfi (2004). 
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In the standard case where polarization phenomena are neglected, the SEE can be obtained 
from Eq. ^ simply by retaining the terms with k - q - Q: 

Jl 

+ YJ^J^f^C,(k^OJ^k^OJ,,T)]p'^(J) (4) 
+ YjCiik = OJ^k = OJi, T)p'^iJ,) + J^Cs(k^OJ^k = OJu, r)p0(7„) 

Jl Ju 

Since, 



C(k^OJ^k = OJ', T) = ^(7 ^ /) 

V2/ + 1 p°(7) = Nj (5) 

where ({ J — > J') is the usual usual fine structure collisional transfer rate and Nj is the population 
of the level J, Eq. Q becomes (with evident notations): 

[j Nj\^^^ = -[ 2 -> D + 2 (,{J ^ Ju, T)] Nj + Y^UJ^ Ju T) Nj, 

Jl J,, Jl 

+ Y,^s(J^Ju,T)Nj, (6) 

J„ 

Therefore, we recover the usual SEE for the level populations. 



3.2. Numerical results 

It is important to note that it is very difficult and sometimes impossible to treat collision pro- 
cesses involving heavy ions like Ba ii by standard quantum chemistry methods, regardless of 
whether they are simple or complex. A semi-classical theory underlying the calculation of the 
depolarization and polarization transfer rates by isotropic collisions of ions with neutral hydro- 
gen is detailed and validated in Derouich et al. (I2004I I. They found that the percentage of error on 
the semi-classical rates with respect to the available quantum chemistry rates is 4% to 7% (see 
Sect. 6 of Derouich et al. 12004b . For computing the collisional rates of Ba ii, we proceed level by 
level since it is not possible to use an unique value of the so-called Unsold energy (Unsold 1927, 
Unsold 1955) as for neutral atoms. Although the fact that this theory is of semi-classical nature, 
the close coupling is taken into account. We apply our collisional numerical code in order to cal- 
culate the collisional transition matrix by solving the Schrodinger equation for the p and cZ-sates 
of Ba II. Then, we obtain the transition probabilities in the tensorial irreducible basis. Afterwards, 
these propabilities are integrated over impact parameters and Maxwellian distribution of relative 
velocities to obtain the depolarization and the transfer rates. We perform calculations varying the 
temperature to obtain the best analytical fit to the collisional rates. 

1 . In the case of the level ^D^ 12, the effect of the collisions is given in the tensorial representation 
by: 



L - (5/2, r)xp|;(7 = 5/2) 
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dpliJ = 5/2) 
dt 



C,{k = 7 = 5/2 ^ ;fc = 7, = 3/2, T) 

Lu = -^^'(5/2, T)xpl(J^ 5/2) 

+[Ciik = 2J = 5/2<^k = 2J, = 3/2, T) 



xp°(7 = 3/2) 



(7) 



(8) 



xp^(7 = 3/2) 



According to our calculations: 

R\5/2, T)is-') = ^Csik = OJ = 5/2^k = 0Ji = 3/2, T) 



-x-x2.88x10-%h(^) exp^i^^-^ 



Ciik = 0J = 5/2^k = 0Ji = 3/2, T)is-^) = 2.88 x lO'^ "h(^) 
R\5/2, T)is-') = 3.44 x 10"' «h(^) +R\5/2, T)is-') 



0.464 



1 \0.26 



Cik = 2J = 5/2^k = 2Ji = 3/2, T)is-') = 1.27 X 10"^ «h(^) 



(9) 
(10) 

(11) 
(12) 



where T and riu are respectively expressed in Kelvins and cm"^ and kg is the Boltzmann 
constant. Given the low degree of radiation anisotropy in the solar atmosphere, one can safely 
solve the SEE neglecting the elements with k > 2, i.e. p^(/ - 5/2). 
2. In the case of the level the effect of the collisions is given in the tensorial representation 
by: 

rdp^iJ = 3/2). 



[- 



dt 



. dpp = 3/2) 
^ dt 



where. 



Icou = -^°(3/2, T)xp°iJ = 3/2) 

+Csik = 0J = 3/2<^k = 0Ju= 5/2, T) x pg(7 = 5/2) 

Icou = -^'(3/2, r)xp2(7 = 3/2) 

+Csik = 2 7 = 3/2 <- A: = 2 7„ = 5/2, T) x p^(7 = 5/2) 



(13) 



(14) 



R\3/2, T)is-^) 



3 , J \0.464 

- x 2.88 X 10-^ nnl 



Csik = 7 = 3/2 <- A; = 7„ = 5/2, T){^~^) 



- X 2.88 X 10" 



.(—1 
V5000/ 

J \0.464 J7 



(15) 
(16) 



exp 



y=5/2 - ii/=3/2 



'""(5000) keT 
R\3/2, T)(s-') = 2.53 X 10"^ «h(^) ' +/?°(3/2, r)(s-') (17) 



Csik = 2J = 3/2<^k = 2J„= 5/2, T)is.-^) 

2 
3 

3. In the case of the level 3/2, one has: 



(18) 



xl.27 X 10-^«h( ^ -^^=5/2-£y=3/2 



''(5000) 



X exp 



knT 



rdpl{J = 3/2). . . 

[-^t L = -^"(3/2, r)xp0(7 = 3/2) 

+[Ci{k = 0J = 3/2<^k = 0Ji= 1/2, 7-)] Xp^iJ = 1/2) 



(19) 
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[^^t ion = -^^3/2, T) xpliJ = 3/2), (20) 

we find that: 

R'\3/2, r)(s-') = ./^ X - X 7.44 X IQ-"^ J—t'' exp ~ (21) 

Ci{k = 07 = 3/2^^ = 07, = 1/2, r)(s~') = 7.44 x 10"^ «h(^) (22) 



/?2(3/2, r)(s-') = 6.82 X 10-%h(^) +^3/2, r)(s-') (23) 

4. In the case of the level the effect of the collisions is: 
rc/p"(7= 1/2), „ 

+[Cs(k = 7 = 1/2 <- ;t = 7„ = 3/2, T)] xp°(7 = 3/2) (24) 

where: 

R\l/2, r)(s-') = y 2 ^•'^^ ""(5000) ^^^^ 

C5(;t = 07=1/2^^ = 07,, = 3/2, r)(s-') = (26) 

1 ,. ,„-9 / 2^ \°'' ____ £y=3/2 - £y=i/2 



(J' \U.:t» 
exp 
5000/ ^ 



2 "\5000/ " /tijr 

5. Finally, elastic collisions with neutral hydrogen do not affect the ground level ^5 1/2]^ 

L dt J™ii ^ ^ 



4. Implication in the Hanle effect diagnostics 

First of all, we mention that, since we are not solving the radiative transfer problem, our results 
have to be considered as complementary informations to the models taking into account radiative 
transfer without collisions and realistic multilevel models. We introduce the collisional effect 
together with radiative rates and we determine the emergent fractional polarization in the limit of 
tangential observation in a plane-parallel atmosphere where the cosine of the heliocentiic angle 

If we calculate the polarization in the framework of the 5 levels-5 lines model as a function 
of the neutral hydrogen density nH (see Fig. [3]), one could remark that colUsions start to play a 
notable role for riu ~ 3 x 10'^ cm"^. The alignment of the metastable levels ^Dj,/2 and '^0^/2 
start to diminish due to collisions and thus the polarization of the Ban line at /14554 decreases. 



^ In linear polarization studies, singlet levels with total angular momentum J= 1/2 are immune to col- 
lisions involving transitions inside the electronic state ~S 1/2- But if the 18% odd isotopes of Ba n are not 
neglected, the levels J = 1/2 are split into hyperfine levels due to coupling with nuclear spin / = 3/2: hyper- 
fine levels F = I and F = 2 can be aligned and can be affected by collisions. 
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Fig. 2. Partial Grotrian diagram of Ba ii showing the levels and the spectral lines of the 5 levels-5 
lines model. Note that the level spacings are not to scale. 

Indeed, for nu ~ 3 x lO'^ cm"^, the ratio of the polarization degree p divided by the zero- 
collisions polarization p,„ax is ~ 0.9 (i.e. a collisional depolarization of ~ 10%). However, in the 
framework of the simplified model (Fig. [T}, a collisional depolarization of ~ 10% of the /14554 
line is attempted only where riu ~ 3 x 10'^ cm"-' (i.e. ~ 100 times larger). This is because the 
upper level ^Pj,/2 of the /14554 line start to be affected by collisions solely for densities Uh > 10'^ 
cm"^ (see Fig.jsjl. 

An estimate height of formation of the A4554 line is h ~ 800 km which corresponds to a 
neutral hydrogen density riu ~ 2 x 10'"* cm"^ and a temperature of the formation of the lines 
T ~ 5350K (e.g. model C of Vernazza et al. 1981). In these typical conditions of formation of 
the Ba n line at /14554, the polarization degree calculated using the simplified atomic model of 
Fig. [T] is practically insensitive to collisions. This may yield, incoiTectly, to the conclusion that 
collisions do not affect the Ba ii line at /14554. In the contrary, we find that, using the 5 levels-5 
lines atomic model, the polarization degree is decreased by ~ 25% because of the collisions. 
Therefore, the collisions are an important ingredient to quantitatively interpret this line. This is 
the main conclusion of the present work. 

The polarization in the absence of collisions is called p,„ax- Since the impact approximation 
is well satisfied for collisions between neutral hydrogen atoms and perturbed ions in the solar 
atmosphere, the collisional rates are simply proportional to the hydrogen density (see their an- 
alytical expressions in Sect. 3). To study the effect of collisions on the polarization degree, we 
report in Fig.[3]the variation of the ratio pipmax as a function of the hydrogen density. In Fig.|3] 
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we make together the results obtained for the simplified model (full line) and for a more realistic 
model (dotted line). We define the percentage of eiTor in evaluating the collisional sensitivity of 
the /14554 line as 



To assess the sensitivity of the microturbulent magnetic field determination to the collisions, 
we proceed as follow: 

1. In the absence of collisions, we determine the polarization degree po which corresponds to 



2. In the absence of collisions, for each value of an eventual collisional inaccuracy Ap deter- 
mined in Eq. (28), the magnetic field B giving the polarization p - po x (l + A/9/IOO) is 
retained. In other words, we introduce a perturbation Ap to the polarization po and we calcu- 
late the corresponding magnetic field. We notice that, in general, the effect of the collisions 
can be overestimated or neglected implying that the polarization should be written generally 
as p = /7oX(l + A/7/100). Here, we investigate only the more typical case where the collisions 
are neglected which corresponds to p = po x (I + A;?/ 100). 

3. The percentage of error in evaluating the magnetic field due to the neglecting of the collisions 
is defined by: 



We obtain AB for each value of Ap (i.e. for each «h)- Afi reaches the asymptotic value 100% 
if = 2Bo or larger. 

The choice of the couple (Bq, po) as a starting point before perturbing the polarization degree 
by a Ap is arbitrary. An other starting couple should give a rather similar AB provided that the 
magnetic field is suitable in the sense that it is well included in the domain of sensitivity of the 
Ban /14554 line to the Hanle eff'ect. 

In Fig. [4j we show the percentage of error AB as a function of «h. For instance, at «h ~ 
2 X 10'"* cm"^, where the neglecting of collisions induces an overestimation of the polarization 
by ~ 25%, the value of the magnetic field is overestimated by ~ 35%. As it could be easily seen 
in Fig. [4] depolarizing collisions can be neglected for nH <~ 10'^ cm"^ since the percentage error 
on the determination of the magnetic field is almost zero. At the asymptotic values of nu (i.e. nu 
>~ 10'^ cm""*), the calculation of the scattering polarization is highly affected by the neglecting 
of the collisions and the information about the magnetic field is completely lost (the percentage 
of eiToris ~ 100%). 

' Bq is the magnetic field strength for which one may expect a sizable change of the scattering polariza- 
tion signal with respect to the unmagnetized reference case. 




(28) 




the critical magnetic field of the Ban /14554 line. Bo - 9 G.|^ 



AB = 1 - (B/Bo) X 100 



(29) 
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Ap n=2^0^^ 

n I 



t 




(cm"-') 



Fig. 3. Emergent linear polarization amplitudes p divided by the zero-collisions polarization p^ax 
versus the density of neutral hydrogen tin- (1) Full line: variation of the linear polarization in the 
framework of the simplified model. (2) Dotted line: same as in (1) but in the framework of the 5 
levels-5 lines model. 
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Fig. 4. Percentage of eiTor on the magnetic field determination due to the neglecting of the colli- 
sions with neutral hydrogen as a function of nn- 



5. Concluding comments 

It is well known that collisions with neutral hydrogen are essential to model the processes gov- 
erning the formation of many polarized lines in the solar photosphere (Di and D2 lines of Nai; 
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Ti I /14536; Sri /14607; etc). We show in this work that, although the neutral hydrogen is about 
ten times less abundant in the low chromosphere than in the photosphere, the polarization of the 
chromospheric /14554 line decreases significantly due to collisions with hydrogen atoms. 

Conclusions about the role of the collisions should not be only based on a simple comparison 
of the inverse lifetime of the upper level of the transition and the coefficients. Although this 
can give useful indications, realistic conclusion should be inferred only from a full introduction 
of the depolarization and colhsional transfer rates in the SEE for multi-level models. 

Lines hke Ban /14554 having upper np-level energy larger than the energy of the (n-l)<i-level 
cannot be treated with a simplified model approximation {n is the principal quantum number). For 
instance, to model quantitatively the K Sr ii /14078 fine one should take into account the ahgnment 
of the J-sates and the infrared triplet Srn /110036, Srn /110327, and Srn /110914. We notice that 
the K Sr n /14078 line was examined by Bianda et al. (1998), but they used the traditional Van der 
Waals approach to calculate colhsional rates and neglected the colhsional depolarization of the 
long hved <i-level. Interestingly, although the Mg n is in the same isoelectronic sequence as Can, 
Sr n, and Ba n, the energy of the upper /7-level of the K Mg n A2196 line is smaller than the energy 
of the <i-level. This is fortunate because the modehng of this fine could be safely performed using 
a simphfied atomic model neglecting the role of the J-level. Furthermore, the K Mg n /12796 line 
is expected to be strongly polarized since the curve showing the anisotropy factor as a function of 
X reaches its maximum around A ~ 2800 A (see Fig. 2 of Manso Sainz & Landi Degl'Irmocenti 
2002). In addition, in given physical conditions, coUisional efl'ects are clearly smaller for Mgn 
than for Srn and Ban. Observations of this fine are difficult from the ground-based telescopes 
but it could be observed with the help of high sensitivity modem polarimeters attached to space 
missions. 

To point out trends in the depth dependence of the magnetic field, Derouich et al. (2006) 
interpreted spatially-resolved observations of the photospheric Sr i /14607 line. In order to extend 
their diagnostic over larger parts of the solar atmosphere, quantitative interpretation of the Ban 
/14554 fine, also observed with spatial resolution, would be highly interesting. To quantitatively 
study this fine, one has to account for partial frequency redistribution efl'ects (e.g. Uitenbroek & 
Bruls 1992 and Rutten & Milkey 1979). Here also colhsions should play an important role since, 
besides their efi'ects on the atomic polarization, they could change the frequency of the Ban 
photons. It remains a challenge to develop a general theory for partial frequency redistribution 
of polarized radiation in the presence of arbitrary magnetic fields and including the effects of 
colhsions in a multilevel picture with/without hyperfine structure. 

Acknowledgements. I would Kke to thank Andres Asensio Ramos and Javier TrujiUo Bueno for provinding 
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improved the presentation of this work. 
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